Optical implementation of Deutsch-Jozsa and Bernstein- Vazirani 
quantum algorithms in eight dimensions 
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We report on a fiber-optics implementation of the Deutsch-Jozsa and Bernstein- Vazirani quantum 
algorithms for 8-point functions. The measured visibility of the 8-path interferometer is about 
97.5%. Potential applications of our setup to quantum communication or cryptographic protocols 
using several qubits are discussed. 



PACS numbers: 03.67.Lx, 42.50.-p 

The last decade has seen the emergence of the field of 
quantum information processing. A particularly promis- 
ing application is the concept of quantum algorithms, 
which allow certain problems such as factorizationjij or 
searching^ to be solved much faster than on a classical 
computer. Another algorithm which we will be inter- 
ested in here is Deutsch's algorithm ||, the first quantum 
algorithm ever discovered, which was later generalized by 
Deutsch and JozsaQ (DJ). The DJ algorithm discrimi- 
nates between a constant or a balanced TV-point binary 
function using one single quantum query, while a classi- 
cal algorithm requires O(N) classical queries. It was later 
on adapted by Bernstein and Vazirani (BV) for efficiently 
querying a quantum database pi 



In the present paper, we report on a fiber-optics im- 
plementation of the DJ algorithm using standard tele- 
com optical components and a single-photon detector. 
The DJ algorithm has already been implemented using 
NMR0 (see also [§ for a NMR implementation of the 
BV algorithm), table-top optics Q (optical demonstra- 
tions of other quantum algorithms also include Grover's 
algorithm |ll|), molecular states (l2), and very re- 
cently using an ion trap(p|. However, our setup sepa- 
rates from these realizations (especially that of ||) on 
several major aspects. First, it relies on guided optics 
components, which makes it unnecessary to perform a 
precise alignment, and it is made robust against phase 
fluctuations by use of an autocompensation technique. 
Second, although it relies on linear optics, our realiza- 
tion is relatively efficient in terms of used optical re- 
sources compared to standard linear optical implemen- 
tations of quantum computation. The central idea of 
such implementations consists in representing the basis 
states of a iV-dimensional Hilbert space by N optical 
paths so that unitary transformations are obtained by 
chaining linear optics components that make these paths 
interfere Such implementations seem, however, 
to be inherently inefficient since the space requirement 
(the number of optical components) and the time re- 
quirement both grow exponentially with the number n 
of qubits (with N — 2") En]. In contrast, in our setup, 



the number of components is kept linear in n, while the 
time needed still grows exponentially. Note that any 
implementation of an algorithm involving an arbitrary 
2 n -point function (also called oracle) does in any case 
require exponential resources to simulate this function. 
Therefore, the linear optical implementation of quantum 
algorithms involving oracles can reasonably be made as 
efficient as any other implementation in this respect. For 
all these reasons, our experimental demonstration works 
with a 8-point (3-qubit) function and might probably 
be extended even further without fundamental difficulty, 
while today's largest size optical demonstrator of the DJ 
algorithm involves a 4-point function^. 

Let us start by recalling the principle of the DJ algo- 
rithm. At the core of the algorithm is the oracle which 
computes a function f(x), where x G {0, 1}™ is an n bit 
string, and / £ {0, 1} is a single bit. The DJ problem is 
to discriminate whether / is a constant or balanced func- 
tion, while querying the oracle as few times as possible. 
A balanced function is such that the number of x's on 
which f{x) — is equal to the number of x's on which 
fix) — 1. Classically, 2 n_1 + 1 queries are necessary in 
the worst case, whereas the DJ algorithm requires a sin- 
gle query as we shall see. In this algorithm, n qubits 
are used, and the basis of the Hilbert space is chosen 
as \x) = \x1X2 ■ ■ ■ x n ) where Xi € {0,1}. The quantum 
oracle carries out the transformation 
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where \y) is an ancilla qubit. By choosing \y) — (|0) — 
\l))/y/2, the action of the oracle simplifies into 



i \ oracle / 

\x) ► (- 



(2) 



since \y) then remains unchanged. 

The DJ algorithm starts with the system in the state 
|0) = 1 00 . . . 0) . Next, a Hadamard transform H is 
applied independently on each of the n qubit. Using 
the definition H\x) = 2~ n / 2 £ 2e {o,i}" (-lf z \z), where 
x ■ z = J^. XiZi mod 2 is the inner product of two n-bit 
strings, we see that the Hadamard transform acting on 



2 



PULSE 


TRIGGER 


GENERATOR 





PATTERN 
GENERATOR 



PHASE PATTERN 

•••>•- 










L.D. 1,55luh 




A 









FARADAY 
MIRROR 



FIG. 1: Fiber optics setup implementing the 8-dimensional Deutsch-Jozsa and Bernstein- Vazirani quantum algorithms. 



the initial state simply yields a uniform superposition of 
all states: H\0) = 2~™/ 2 J2 X \%)- This state is then sent 
through the oracle whereupon it becomes 

2 -n/2J2(-l) f W\x). (3) 

x 

The superposition principle allows the oracle to be 
queried on all input values in parallel. A second Hada- 
mard transform is then carried out to obtain the state 

2- n Y^(- l Y' z+f(x) \ z ) ( 4 ) 

which is finally measured in the z basis. One easily de- 
duces from Eq. (|4|) that when / is constant, the prob- 
ability of measuring |0) is one. In contrast, when / is 
balanced, this probability is always zero, so the DJ algo- 
rithm can distinguish with certainty between these two 
classes of functions by querying the oracle a single time. 
The BV variant of this algorithm is also based on the 
transformation leading to Eq. (^). Suppose that the or- 
acle is restricted to be of the form fj (x) — x ■ j where 
j € {0, 1}" is an arbitrary n-bit string. The aim is to find 
the bit string j with as few queries as possible. Classically 
one needs at least n queries since each query provides one 
independent bit of information at most about j. Quan- 
tum mechanically, a single query suffices since using fj (x) 
in Eq. (^) shows that the measurement outcome is z = j 
with probability. 

Our all optical fiber (standard SMF-28) setup is il- 
lustrated in Fig. 1. Initially, a 3 ns light pulse is pro- 
duced by a laser diode at 1, 55 /im, attenuated by an op- 
tical attenuator (Agilent 8156A), and then is processed 
through three unbalanced Mach-Zehnder (M-Z) interfer- 
ometers with path length differences A; (I = 1,2,3) obey- 
ing A3 = 2A2 = 4Ai = 15 ns. Each MZ interferometer 
doubles the number of pulses so that, at the coupler C6, 
we get eight equally spaced pulses. This corresponds to 
the action of the Hadamard transform on the three input 
qubits in state |0). The pulses are then reflected by a 
Faraday mirror and, on their way back, are modulated 
by a phase modulator (Trilink) commanded by a pat- 
tern generator (Agilent 81110A) which selectively puts 
a phase shift of or ir on each pulse according to the 
8-point function [see Eq. (||)]. The pulses then pass back 



through the three M-Z interferometers, thereby realiz- 
ing a second triple Hadamard transform, and are sent 
via a circulator to a single-photon detector (id Quan- 
tique id200), which completes the DJ or BV algorithm. 
The additional delay lines L\, L2, and L3, which obey 
L3 > 2L2 > 4Li > 8A3, ensure that the different out- 
puts of the DJ or BV algorithm all reach the detector at 
different times. The delay lines L2 and L3 also contain 
an isolator so that the pulses are not transmitted on their 
way to the mirror. The photodetector was gated during 
5 ns around the arrival time of each pulse by the pattern 
generator. The output of the detector was registered us- 
ing a time to digital delay converter (ACAM-GP1) con- 
nected to a computer. All delays Aj and Lj were chosen 
to be integer multiples of Ai within 0,2 ns. All electronic 
components were triggered by a pulse generator (Stand- 
ford Research Inc. DG535). In order to maximize the 
visibility, polarization controllers were introduced in the 
long arm of each M-Z interferometer and in front of the 
polarization-sensitive phase modulator. Once optimized, 
the setup was stable for days. 

This implementation of the DJ and BV algorithms dif- 
fers from an earlier optical implementation of the DJ 
algorithmQ in several important aspects. First, we run 
the algorithm for n = 3 qubits and, more importantly, 
we measure all 8 outcomes, which makes it possible to 
realize the BV algorithm as well (the previous implemen- 
tation |j| works with n — 2 qubits and only measures the 
outcome z — 0). Second, we operate at telecom wave- 
lengths in optical fibers using a setup closely inspired 
from the "plug-and-play" quantum cryptographic sys- 
tem developed by Gisin and collaborators (see e.g. p7[). 
For this reason, the present setup in a slightly modified 
version can be adapted to implement quantum cryptog- 
raphy using higher dimensional systems or to illus- 
trate quantum communication complexity protocols |fl9| 
over distances of a few kilometers. These potential ap- 
plications will be discussed below. Third, the used re- 
sources are quite different in the two implementations: 
when scaled to a large n, the implementation of |j) re- 
quires exponential time and exponentially many optical 
elements, whereas our implementation also requires expo- 
nential time but only a linear number of optical elements. 
This is because the n qubits are realized as 2™ separate 
optical paths in |l, whereas in our case they are repre- 
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sented as 2 n light pulses traveling in a single optical fiber, 
extending naturally the "time-bin" realization of qubits 
used in fTij] . The small number of optical elements in 
our setup therefore implies that it is relatively easy to 
increase the number of qubits n while keeping the opti- 
cal setup stable. As we shall see, a disadvantage of our 
setup is that the Hadamard transform can only be im- 
plemented with a probability of success of 1/2. Since 2n 
Hadamard transforms are needed for the DJ algorithm 
with n qubits, the resulting attenuation is 2~ 2n . 

Let us now prove that our optical setup indeed realizes 
the DJ and BV algorithm. The quantum state describing 
the eight pulses at coupler C6 can be written as 



111 

oc exp 

x=000 



i 2^{kxi A; + ttxi) 



1=1 



where x stands for {x\,x>},Xz) and \p) denotes a pulse 
located at position p. The three bits x±,X2,X3 = 0,1 
label whether the pulse took the short (x = 0) or the 
long (x — 1) path through each interferometer. The fac- 
tor exp[ifc^ ; xiAi] with k being the wave number takes 
into account the phase difference between a pulse travel- 
ing along the short or long paths of the interferometers. 
The factor exp[z ^ nxi] takes into account the phase ac- 
cumulated at the couplers: if the pulse takes the short 
path, it is transmitted at two couplers, whereas, if it 
takes the long path it is reflected twice. After reflection 
at the Faraday mirror and phase modulation, the pulses 
cross again the three M-Z interferometers and reach the 
photodetector in the state 
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| +w)Aj+*Li)) (6) 

i=i 

where the bits 2/1,2/2,2/3 equal or 1 according to whether the pulse passed through the short or the long path of 
each of the interferometers on its way back, and the bits z±, z%, 2:3 equal or 1 according to whether or not the pulse 
exited each of the interferometer in the path containing the delay lines L\, L 2 or L3. Note that we put Z4 = 0. We 
have again taken into account the phases induced by transmission or reflexion at the couplers C1-C6. The final state 
contains 120 pulses, but we are only interested in the eight pulses such that X\ + y\ — x-2 + 2/2 = ^3 + 2/3 = 1, which 
are those that exhibit 8-path interference. The other pulses are filtered out in the computer analysis (they correspond 
to different time bins). The final state then becomes 

|V) OC ^H) Z1 (-1) Z2+Z3 ^(_l)/(*i.»a,«s)+«i(*i+«)+» a (* a +«8)+«a*a | Al + A 2 + A3 + Z X L X + Z 2 L 2 + Z 3 L 3 ) (7) 
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( K x i + vi)&i + n ( x i +vi- yi( z i + z i+i) h ^— — 



By relabeling the time bins according to the substitution 
z\ — > z\ + Z2 + Z3 mod 2, Z2 — ► z 2 + Z3 mod 2, and 
z 3 — > Z3, this equation coincides (up to irrelevant phases 
and an overall normalization factor) with Eq. (|) with 
the 8 logical states \z) identified as specific time bins. 

The setup thus realizes the DJ or BV algorithm, the 
main difference with an ideal algorithm being an extra 
attenuation by a factor of 2 -7 . A factor 2 -3 originates 
from the couplers C2, C4, and C6 because each time a 
pulse passes through these couplers it only has a prob- 
ability 1/2 of exiting by the right path. Otherwise, it 
is absorbed by the isolators II and 12 or by the uncon- 
nected fiber pigtail at coupler C6. Another factor 2~ 3 is 
due to the filtering out of the 112 pulses produced on the 
way back that do not correspond to 8-path interferences. 
The remaining factor 2~ 1 is due to the coupler C7. This 
overall loss of 21 dB could be remedied by replacing the 
couplers C2, C4, C6, and C7 by optical switches which 
would direct the light pulses along the appropriate path. 
High speed, low loss optical switches are not available 
commercially at present, so we had to use couplers in the 
present experiment. However, we emphasize that this is 



a technological rather than a fundamental limitation. 

In order to characterize the performances of our setup, 
we considered the 2" oracles of the form fj(x) — x ■ j 
and fj(x) = x ■ j + 1 mod 2 (i.e., the oracles in the BV 
algorithm and their complements). For each oracle fj (or 
/„•), we ran the algorithm 500,000 times and registered 
the number of counts in time bin z, denoted as Nj(z) (or 
Nj(z)). The algorithm gives constructive interference in 
the time bin z = j for the oracle fj or fj , and destructive 
interference elsewhere. We then computed 

A) 2{n z (z) + N j (z) + N z (z) + N 3 (z)J {> 

for each pair of oracles with j ^ z, and calculated the vis- 
ibility V(z) in time bin z by taking the average of Vj(z) 
over all values of j. The measured visibilities V(z) for 
2 and 3 qubits are shown in Table I. Remarkably, they 
remain relatively high when going from 2 to 3 qubits in 
spite of the fact that 8 path interferences are involved. 
This is because the path differences are automatically 
compensated and only n + 1 polarizations must be ad- 
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_z 1 2 3 4 5 6 7 8 

V(z) n = 2 98.4 97.4 98.5 98.6 

V{z) n = 3 96.78 97.99 97.68 97.32 97.33 97.56 97.37 97.28 

TABLE I: Measured average visibility V(z) in the zih time 
bin for the DJ or BV algorithm with n — 2 and n — 3 qubits. 

justed. It should therefore be relatively easy to go beyond 
n = 3 without significantly decreasing the visibilities. 

Because of the attenuation in our setup along with the 
quantum efficiency (~ 10.5 %) and the dark count prob- 
ability (~ 10~ 4 ns _1 ) of our detector, the signal-to-noise 
ratio was not high enough to perform these visibility mea- 
surements in the single-photon regime. For this reason, 
in our experiment, approximatively 20 photons per pulse 
entered the phase modulator (oracle) in 4 dimensions, 
and approximatively 50 in 8 dimensions. However, mini- 
mal modifications should allow us to decrease the number 
of photons while keeping the signal-to-noise ratio con- 
stant. In particular, by reducing the pulse length or us- 
ing two detectors instead of the coupler C7, it should be 
possible to operate in the single-photon regime in 4 (and 
possibly 8) dimensions. Moreover, as already mentioned, 
the Hadamard transform could be rendered deterministic 
by using fast low-loss optical switches instead of couplers, 
which would strongly reduce the losses. 

The present experiment can be extended in several 
ways. For example, one could implement the distributed 
Deutsch-Jozsa problem Jl9[ where two parties, which re- 
ceive each as input a 2"-bit string (denoted as / and 
g), must decide whether / = g or / differs from g in 
exactly 2™ _1 bits (they are promised that only one of 
these two cases can occur). The value of n for which the 
gap between the classical and quantum algorithms sets 
in is unknown, but recent results suggest that it could 
be as early as n = 4 J2(J. The distributed Deutsch-Jozsa 
algorithm could be realized with a slight modification of 
our setup in which the two parties, separated by by a few 
kilometers of optical fiber, would use each a phase modu- 



lator. This quantum protocol can also be easily adapted 
to multidimensional quantum cryptography, where two 
parties randomly choose their patterns / and g. By pub- 
licly revealing part of / and g, the parties can use the 
correlations between the measurement outcomes to es- 
tablish a secret key. A detailed analysis shows that for n 
= 2, this exactly coincides with the 4-dimensional cryp- 
tosystem based on 2 mutually conjugate bases, which has 
been shown to present advantages over quantum cryptog- 
raphy in two dimensions pi]] . A final potential applica- 
tion of this setup is to test quantum non locality using 
the entanglement-based Deutsch-Jozsa correlations p2| . 
An entangled state of 2™ time bins must be produced, 
for example using the source p3| , and each party must 
then carry out phase modulation and a Hadamard trans- 
form (which we have demonstrated are easy to realize on 
time bin entangled photons). The correlations between 
the chosen phases and the measured time of arrival of 
the photons at each side should exhibit quantum non lo- 
cality. It has recently been shown that these correlations 
are non local for n — 4 pb]| , and that they exhibit ex- 
ponentially strong resistance to detector inefficiency for 
large n ]24| , which means here that the 3 dB losses at 
each Hadamard transform can in principle be tolerated. 

In summary, we have demonstrated, by implementing 
the DJ and BV algorithms for 3 qubits, a simple and 
robust method for manipulating multidimensional quan- 
tum information encoded in time bins in optical fibers. 
We anticipate that our method will have wide applica- 
bility for quantum information processing and quantum 
communication using higher dimensional systems. 
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